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Abstract 

Massless and massive scalar fields and massless spinor fields are considered at arbitrary 
temperatures in four dimensional ultrastatic curved spacetime. Scalar models under con- 
sideration can be either conformal or nonconformal and include selfinter action. The one-loop 
nonlocal effective action at finite temperature and free energy for these quantum fields are 
found up to the second order in background field strengths using the covariant perturbation 
theory. The resulting expressions are free of infrared divergences. Spectral representations 
for nonlocal terms of high temperature expansions are obtained. 



PACS: 04.62,+v, ll.10.Lm, ll.10.Wx 



1 Introduction 



Finite temperature field theory has been developed in a series of seminal papers 
[|, [|. Nowadays it is an actively growing branch of theoretical physics 0. Ther- 
modynamical properties of thermal quantum fields in the presence of background 
fields are very important for a large number of applications in high energy physics, 
astrophysics, and cosmology. However most of these studies are devoted to the situa- 
tion when background fields are constant (homogeneous) 0, . This particular form 
of the effective action, the effective potential |7|, §, when large background fields are 
taken into account nonperturbatively, is useful for study of phase transitions in the 
early Universe or quark-gluon plasma. For a long time, the opposite situation, when 
background fields are small but rapidly fluctuating, lacked investigation even in zero 
temperature field theory. Traditional tools of quantum field theory, like the short 



proper time Schwinger-DeWitt expansion |9], [T(| pT| , are intrinsically local, hence, 
they miss nonlocal contributions. As a consequence of this deficiency artificial in- 
frared divergences appear in the perturbative effective action for massless fields, and 
perturbation theory breaks down. Finite temperature effects also contribute to in- 
frared divergences [Q, an d methods of diagram summations have been developed to 



improve the perturbation series |L2], pL3| . 

To deal with massless field theories properly, such as gauge field theories or quan- 
tum gravity, Vilkovisky suggested a new powerful method [TJ] which is known as the 
covariant perturbation theory |15[ [16, [17|, [lSfl . In these papers it was shown that 
infrared divergences are artificial and brought into existence by a mode of calculation 
rather than by a field theory. They disappear after summation of terms with infinite 
number of derivatives acting on background fields, which results in nonlocal terms 



entering the effective action [fL5 |. Such a summation can only be performed in a given 
order in background field strengths. 

Thermodynamics of an ensemble of quantum fields in equilibrium in static curved 
spacetimes is well defined, and most properties of the system can be derived from 
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its free energy ||19| , |20| , |2jJ . In this paper we consider ensembles of scalar and spinor 
fields in the presence of external ultrastatic gravitational field. The scalar models 
may have an arbitrary interaction potential and an arbitrary coupling to gravity. We 
employ the method of covariant perturbation theory to find the finite temperature 
effective action and the corresponding free energy of these quantum fields on highly 
inhomogeneous gravitational backgrounds. An example of the situation when finite 
temperature effects on curved background are important, and, thus, nonlocal effective 



action is needed, is the Hawking radiation by black holes [£2], p23 |. 

The paper is organized as follows. In the next section we describe how to obtain 
nonlocal free energy at finite temperature with help of the covariant perturbation 
theory. In section |3| we derive the free energy of interacting massless scalar fields 
and study its high temperature behavior. Massive scalar fields at high temperatures 
are treated in section (|. In section |5] we derive free energy for massless spinor fields 
at finite and high temperatures. Conclusion and discussion of possible applications 
and extensions of obtained results can be found in section ||. We place necessary 
complicated computations into appendices [A] and [B]. 

2 One-loop effective action and free energy of quan- 
tum fields in ultrastatic spacetimes 

Let us consider fields (p described by the classical action S(<p) and the corresponding 
canonical Hamiltonian in a generic curved static spacetime. Statistical free energy 
Fs of the ensemble is defined as the trace of logarithm of eigenvalues of the normal- 
ordered Hamiltonian. In canonical quantization scheme, ultraviolet divergencies are 
traditionally subtracted from Fs by the normal ordering prescription plj [24]. On 
the other hand, in the imaginary time formalism of Matsubara the problem of 
finding free energy of the system in equilibrium reduces to the computation of the 
path integral of the Euclidean quantum field theory and the corresponding effective 
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action W p [g] 



The temperature T = 1/(/cb/3) enters the calculation via the condition of (anti)periodicity 
in the Euclidean time r imposed on quantum fields with (fermi) bose statistics 

<p(x,T) = ±<p{x,T + 0) (2) 

(the Botzmann's constant fee = 1 everywhere). 

The canonical free energy F$ and the thermal renormalized Euclidean effective 
action W@ are closely related to each other and differ [TP], only by terms F that 
are independent of temperature 

iWM = *SM + FM, (3) 

where <p —< (p > are mean fields. The effective action is usually regularized us- 



ing covariant methods, e.g., zeta-function [25|, |26[, dimensional []T0| , etc., while the 
canonical free energy is regularized via normal ordering of operators. The difference 
F[(j)\ is related to different ways of taking into account vacuum energy contributions 
in covariant and canonical regularization schemes. The covariant approach is more 
appropriate to our problems since it is consistent with calculations of the stress tensor 
and vacuum polarization effects in external fields. In any case, it is easy to compute 



F((f)) which is temperature independent and local f24 |. Henceforth, we restrict our 
consideration to calculation of the one-loop Euclidean effective action W 13 and the 
corresponding covariant Euclidean free energy 

Ff> EE IW 13 . (4) 

P 

We calculate the free energy of quantum fields on static background fields which 
include mean field (f> and static gravitational field. The Tolman temperature of such 
a field system in equilibrium is not constant throughout the static space. It is more 
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convenient to perform calculations of temperature effects in the Euclidean ultrastatic 
(optical) spacetimes, 

ds 2 = g^dx^dx" = dr 2 + g ij dx i dx j , (5) 

where local temperature is constant throughout the space. Ultrastatic and static 
spacetimes are related to each other by a conformal transformation of the metric. 
Conformal properties of the effective action have been studied in detail p7[ , and 
applied to free energy calculations by Dowker and Schofield [25]. Using scaling 



properties of the finite temperature zeta functions it was shown that the difference of 
free energies in two conformally related spaces does not depend on temperature. Then, 
all temperature dependent terms can be found from the free energy in an ultrastatic 
space, where solution of the problem simplifies significantly. This difference can be 
found by integrating the conformal anomaly, but the method of Ref. [2£| works for 
generic nonconformal operators as well. 

The brief outline of our research program is to calculate free energy in an ultra- 
static space, and then using a relation between free energies in static and ultrastatic 
spaces to express the final result in terms of quantities defined in a physical (static) 
spacetime. In this paper we implement the first and most complicated step of obtain- 
ing W 13 and F@ on the ultrastatic metric (|5|). 

Let us consider quantum n-component scalar field <p = <pa, i = 1 . . . n, which 
satisfies the equation, 

[(□-ifl)i + P(0]¥> = O, (6) 



Our notations correspond to those of Refs. [10], [IJ|: the Laplacian □ = g >JjU 'V ll 'V u 
is constructed of covariant derivatives which are characterized by the commutator 
curvature 

(V M V, - V,V M V = (7) 

This quantity is, of course, zero for scalar fields, but we will need it in section ||] where 
spinor fields are considered. The potential P may depend on the metric and mean 



field 0, which is a part of classical background. Thus, this class of models includes 
self-interacting fields. The vanishing potential P = corresponds to the case of free 
conformal scalar fields, and P = 1R/6 - to the minimally coupled free scalar fields. 
The overhat symbol indicates the matrix structures, P = P A b, and term R/6 in (^) 
is explicitly singled out for convenience. The three field strengths R^ u , TZ^, P will 
be also referred to as curvatures. This massless field theory will be generalized to the 
case of massive fields in section f|. 

The one-loop Euclidean effective actions W@ is defined in terms of the functional 
trace of the heat kernel, 



Ir^TrK^s), (8) 
2 Jo s 



where the heat kernel K^{s) is the periodic in Euclidean time solution of the problem 

1 



d 



l(U - -R(x)) + P(x) \K p (8\x,xf) = l5(8)6(x,x f ), (9) 
6 J 



') = K^{s\t + P, x; r', x'). (10) 
The functional trace is understood as 

TrK(s) = [ d D xtr K(s), (11) 



with the tr standing for the matrix trace, e.g., trl = 5 a a, trP = P a a- 

The thermal (periodic in the Euclidean time) heat kernel K@ can be expressed as 



an infinite sum of zero temperature (vacuum) heat kernels H, |30] 



K p {s\t,x;t',x') = £ K{s\t,x;t' + (3n,x) (12) 

n=— oo 

This image sum is equivalent to summation over Matsubara frequencies in a momen- 
tum space representation in thermal field theory. The image sum in the context of 



Casimir energy calculations was introduced in Ref. |3T| . 

Temperature effects are inherently connected to the imaginary time. It is conve- 
nient to factorize the heat kernel into temporal and spatial parts, 

k(s\r, x; r>, x') = jJ^p exp (-^l)k^(s\x; x% (13) 
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which is possible to do in ultrastatic spacetimes. Then, the trace of the heat kernel 
takes a form [ [TT| , 

TtK^s) = e 3 (Q,e-£)j^—Jd 3 xtrK®(s\x,x), (14) 



when expressed in terms of the Jacobi theta function ||32|| , which is defined in a usual 

way, 

n=oo 

9 3 (a,b) = £ e 2nai b n \ (15) 

n=— oo 

The free energy of quantum fields in static spacetime F@ is defined via the finite 
temperature Euclidean effective action W 13 and can be written in the form, 

**= -±r-TrK0(8). (16) 
zp Jo s 



The vacuum mode n = in the infinite sum ([14])-(|15D corresponds to the zero 
temperature effective action which suffers ultraviolet divergencies [IT], Fortu- 
nately, this is the only divergent term of the sum so it is convenient to treat it 
separately We subtract the zero temperature (/3 = oo) free energy F°° from F 13 and 
renormalize it with the use of the zeta function regularization [SSL ISJ: 



2de 



fi 2t [°° ds 

Jo —^ lK(S \ e = 

where /x is a mass-like regularization parameter and T is the gamma function. Ff° 
will be combined with n ^ terms at the end of our derivations. Therefore, we 
compute 

FL-FZ= "2/0 -(^e-^-^p^/d^tr^ 3 )^;^). (18) 

The heat kernel (s) is defined as a solution of Eq. (§) with the three dimensional 
operator, 



He) 



(17) 



1C11 



d_ 
d7 



1 



1(A - + P(x)\ I K {3) (s\x; x) = I5(s)8(x, x'). (19) 



In this case the three dimensional Laplacian A, and potential P(x) and the curvature 
Rij{x) are defined on a three dimensional hypersurface r = const of the ultrastatic 
spacetime. 



Many methods have been developed for calculation of the trace of the heat kernel 
0, ^4]. Most of them (see reviews [11, |HJ) reduce to various representations of its 
small s expansion, 



Trir( S ) = Ii ^/dxV /2 (x)E^r 



0. 



(20) 



n=0 



However, as soon as the inverse temperature (3 is finite, the behavior of the heat kernel 



at large values of proper time s becomes very important fll9[ . Therefore, expansion 
( p0|) is not suitable for our task of finding the free energy at finite temperature. 
Besides, Schwinger-DeWitt coefficients a n are local functions of background fields, 
henceforth, nonlocal free energy cannot be derived using (|20|) . To solve the problem of 
obtaining nonlocal free energy at finite temperature we have to resort to the covariant 



perturbation theory [TJ|, [L6|, [T7|, [T^]. There is no need to repeat derivations of the 
covariant perturbation theory here because an expression for Tr K(s) is already known 
in arbitrary D-dimensions [[16|, fjl| . In this paper we will take it up to terms quadratic 
in curvatures, 



R^h{-sU)R^\ + Rf 2 {-sU)Ri + Pf 3 {-sn)R 



+ Pf 4 {-sU)P + UM-8U)W a '] } + 0[3? 3 ]. 



(21) 



Analytic functions fa (form factors) have the dimensionless argument sD. (The ap- 
pearance of nonlocal form factors in the momentum space representation of the effec- 



tive action originates in the classical paper of Schwinger fl35|). The form factors act 
on tensor invariants constructed of the set of field strengths R al3 , P, TZ^ character- 
izing background. The collective notation 3? will be used for these curvatures. First 



two terms of the sum (21) are purely local and coincide with first two coefficients of 
the short proper time expansion (p0| ). Formally, the expansion (21) is valid only in 
an asymptotically flat Euclidean spacetime with the topology R D . All background 
curvatures 5i are supposed to vanish at spacetime infinity fl6| . Since we use a pertur- 
bation theory, all calculations are carried out with accuracy 0[3? n ], i.e., up to terms 



of n-th and higher power in the curvatures 3?. The very structure of this curvature 
expansion restricts its validity to background fields satisfying the relation, 

VV3? » K 2 . (22) 

Physically it means that gravitational fields are small in magnitude but quickly os- 
cillate. 

All form factors in ( pT|) can be expressed in terms of one basic form factor 

f(-sD) = ['da e a(1 - a)sn . (23) 
Jo 



Their explicit form reads |T(J 



/I( _ SQ) . MzizifE, (24) 



M-sn) 



l 



36 M J 3 sD (sD) 2 



(25) 



/!(-«□) = ^ /(-*□), (27) 

«--) - l^i- W 

Even though, in the following consideration general covariance is broken because 
of the presence of temperature, we will refer to this curvature expansion as to the 
covariant perturbation theory. In spatial dimensions the covariance remains explicit. 

A few words about validity of this approximation are in order. Since we consider 
quantum fields at some fixed temperature, one can say that the field system in ques- 
tion represents a canonical ensemble. To define a canonical ensemble rigorously we 
have to assume that the fields are in some cavity of a finite volume, as it is usually 



assumed in the presence of a black hole [21|. This assumption should be reconciled 



however, with our method of computation described above, which in the present form 
works only for asymptotically flat spacetimes and requires vanishing background fields 
at spacetime infinity. It is important to note that background field strengths, sources 
of vacuum polarization, have a compact support on a manifold, thus, providing an 



effective volume cut-off. In regard to gravitational field this property is due to the 
presence of the Ricci tensor rather than the Riemann tensor [16]. 



3 Free energy of massless scalar fields 



Let us now compute free energy fll%p of massless scalar fields at finite temperature. 



This case was briefly reported in our letter [36]. After introducing a new variable 
y = (3 2 /4s, first two terms of the trace of the heat kernel (|21|) take the form of the 
integral, 



l o °° dy [03(o,e-») - l] y a - 1 = 2((2a)T(a) 



(29) 



where ( is the Riemann zeta function, T is the gamma function. When a is taking 
values 2 and 1, expression (JE^) gives for the zeroth and first curvature orders coeffi- 
cients 7r 4 /45 and 7r 2 /3 correspondingly. These local contributions to the free energy 
are well known [|19[ |37j and coincide with the first two terms of high temperature ex- 
pansion. Since all information about temperature is separated from tensor invariants, 
we can write down an anticipated form of free energy up to second order in the field 
strengths, 

1 



poo 
ren ren 



+ 



d 6 xg 
1 



^trj 



7T 



-P 



90/3 4 2A(3 2 ' 
Rijli{-A)R ij + R^(-A)R 



327T 2 

+ P 1 S{-A)R + P 7 4 (-A)P] + 0[K 3 ]}. 



(30) 



Then, the problem with the second curvature order is reduced now to calculation 
of the thermal form factors, 



r°° ds 






3 (o,e-£) -1 


/o s 





fi(sA), 



(31) 



where are given by (p^)-(|27|). We show how to compute ([H]) when fi(-sA) is the 
basic form factor (^). After substituting fl2"3] ) into ( PT[ ) and writing down the theta 
function (|l~5| ) explicitly we get, 



7(/?V = A) = 2 £ ['da T^exp 



yn 2 -±-a(l-a)(3 2 (-A)). 



(32) 
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Integration over y produces the modified Bessel function of the second kind 

oo r i 



7(2;) — 4 / daK (nz\J a(l — a) 



(33) 



where z = f3^—A. Change of variables, x = 2ya(l — a), allows us to express 
in terms of the exponential integrals, 

1 



dx 



x 



--Kn 



,nz x. 



x A 



nz 



EU — ) e - nz / 2 -Ei(-—)e nz / 2 



(34) 



Now we can use for the right hand side of (|34T ) its standard form in terms of elementary 



functions [321 and obtain 



roc x 

j(z) = 4 / dt £ 

J0 n=l 



sin(i) 



t 2 + n 2 z 2 /A' 

The sum over n can be evaluated p2| , 

~2vr 1 



(35) 



y i = i 

^ 1 t 2 + n 2 z 2 /A 2 



1 



zt th(2wt/z) t 2 



(36) 



and the resulting expression reads, 



/■OO 

7(2) = 2 / dt sin(t) 
Jo 



2tt 



(37) 



th(2vrt/z) t 2 

As can be seen from (|24])-(pTD, there are two other types of basic thermal form 
factors (with one and with two subtractions). Their derivations can be found in 
appendix A. Applying results (p7[), (|105|) , and ( |1 14|) to the table of form factors we 
obtain for all thermal form factors the following expression, 



7i(/?v^A) = [°°dt 9l {t) 



2tt 



Py/=Et th(27rt/(/3v^A)) t 2 _ 
and gi (i — 1, ...4) are simple combinations of elementary functions 



(3* 



9i(t) 
92(t) 

9i(t) 



2 

1/1 . . COSH j 

— (-sin (t +2 — + 

48^3 W t 



1 ^sin(t) rcos(t) sin(t) 
~\~ 3 



f 1 



t 3 



cos(t) sin(t) i-cos(t) sin(t) 

~T + ~C 2 ~ + \-~ W 



1/1 . . cos(t) sin(i) 

-f-sn(t) + — — — y- 

2V3 KJ t t 2 



2^3 
sin(t) 



(39) 

(40) 

(41) 
(42) 
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The final result for renormalized free energy at finite temperature Ff en is presented 
by a sum of Eqs. (|30|)-(f4T|) and renormalized free energy at zero temperature F l 
After the zeta regularization flT7|) , the latter one takes the form 
1 



32vr- 



J ^xg 1 'hi:{R ijll {-A)R ii + Rj 2 (-A)R 



+ P l3 (-A)R + P 74 (-A)P + 0[K 3 ]}. 
where zero temperature form factors 7»(— A), i — 1, ...4, are 
7i(-A) 



1 

60 



72 (-A) = — 



15 

, Ax 97 



Ts(-A) 
74(-A: 



-M-4U2 



2/ 



(43) 

(44) 
(45) 
(46) 
(47) 



This expression differs from the one obtained using dimensional regularization only 



by unessential additive constants [16|. 

Formulae (|30"D-(fi7]) , we have obtained, are valid at arbitrary finite temperature. 
Now we would like to study asymptotic behavior of the free energy ([Hp in high tem- 
perature regime, the most interesting and the best studied limit. In the framework of 
perturbation theory, the problem boils down to finding j3 — > asymptotic of thermal 
form factors (PB[). We have to be careful while dealing with mutually compensating 



singularities. After relatively straightforward calculations the outcome for (|37| ) is, 

2vr 2 



t(/V-a) 



+ 2 



C(5) 



0. 



(4f 



where C is Euler's constant and ( is the Riemann zeta function. 

Now, expressions for the vacuum free energy fl4~3"D and the high temperature ex- 



pansion of (|30|) match, and can be combined into a single formula. The resulting 
T — > oo expansion of the renormalized one loop free energy takes a form, 

FL = - [d'xg^tr 



7T 2 . 1 * 

1 + 



90/5 



24/3 2 



12 



+ 



1 



lo 1 



32/5 

2l ln (— 



R 



i.i 



25 



-i? 



1 



-.R 



P- 



1152 v^A 12 V^A 



1 =* + P 1 



P 



16tt 

+ /3 



47T 



c 



60 J 



—Pi? + -PP 
180 2 



2 C(3) r 1 



+ (3 



1287T 4 

4 3C(5) r 1 



840 3 



—^—RAR + —PAR + —PAP 
3780 180 12 



R A 2 R ij 
15120 13 



1024tt 6 
+ 0[K 3 ]+0[/3 6 ]}, p 



-^—PA 2 R + —PA 2 P 
1260 + 120 



0. 



(49) 



All local terms of this result perfectly reproduce those of Refs. [|19|, ^8 . The combina- 
tion of quadratic in curvatures terms at the logarithm is just the trace of the second 
Schwinger-DeWitt coefficient a 2 , taken with Riemann curvature expressed via Ricci 
one p, |H], [16j . Higher powers of (3 in Eq. (|49D are also quadratic in curvatures parts 
of 03 and a 4 Schwinger-DeWitt coefficients ^] . 

We obtained the explicit form of all nonlocal terms of the second order in cur- 
vatures. They are proportional to 1//3, and were known to exist p8| . The general 
structure of nonlocal terms is !R I — and, therefore, techniques based on local 



v A 

(small s) expansions could not generate them. Terms of higher orders in curvatures 
fl"8| , |38| will also give nonlocal contribution linear in temperature. 

The meaning of nonlocal structures can be understood from spectral representa- 
tions in terms of massive Green functions |39], pj], [TSfl . For this particular form we 
have the following spectral formula, 
1 2 r°° , 1 



V=A 



dm 
ir Jo m 



2 - A' 



(50) 



A remarkable property of the expression ( |49D is that it contains the only kind of 
nonlocality, (|50|). All logarithmic nonlocalities ln(— A), that are present in F™ n and 
Ff en , have mutually canceled, leaving logarithm temperature dependance in the form 
of ln(/3/i), This local combination is well known in both flat ||, 0] and curved [19 



space thermal field theory. The In (—A) disappearance is still being analyzed in a 
different physical language and in a different setting [jlT] . 

Of course, we are not completely satisfied with the integral representation for 
the free energy at finite temperature (|30|) . Although, it admits a closed form, we 
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would prefer to see F@ expressed entirely in terms of analytical and special functions. 
Indeed, it is possible to obtain such a form after applying the Poisson resummation 



32, 



E e-fc" = ^ E • (51) 

P k=-oo 



n=— oo 



Then, the following identity holds, 




# 3 (0,e 4.")-! = — — E e ^ - /_ dKe ^ ■ (52) 



We compute now the basic thermal form factor, ([31] ) with (|23|) , using this identity 
and separating k = term out of the sum, 

7(/? v- A) = w^ + pV^ L?! arctan ^^ " dKarctan ii^r) ) - (53) 

The fc = mode of the sum gives precisely the leading infinite temperature contribu- 
tion, while the rest can be calculated by employing the following sum, 

g(— (s)-!)-Hi^H 

Adding up the regularized zero temperature form factor, 

7 (-A) = -ln(^)+2, (55) 
we obtain an expression which is valid at any temperature, 

^^ = ^7TK + 1T7TK^ ~ + 21n S- ( 56 ) 



Vr(i-i^' 



Besides an obvious advantage of Eq. fl56|), namely, that it is the formula in terms of 
usual elementary and special functions, the leading infinite temperature contributions 
are present here explicitly. Taking (3 — > oo and (3 —>■ limits, one can readily find 
zero temperature fl55]) and high temperature (|48|) asymptotics of this basic thermal 
form factor. In fact, one can see that the logarithm of the gamma functions' ratio in 



the main result (56) is a sum of all positive powers of (3 in the high temperature limit 
([48|). Hence, it gives a partial (in the given curvature order) summation formula for 
the j3 — > series |[19|| . Eventually, one has to transform fl56|) into a spectral form, the 
procedure we can complete again only at high temperatures, ([501) . This is the reason 
why we refrain from deriving the total free energy in this new representation. 
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4 Free energy of massive scalar fields 

The use of the curvature expansion is crucial for derivation of the massless field 
free energy, because it allows one to avoid artificial infrared divergences. Two other 
advantages of perturbation theory, namely, that free energy can be found at arbitrary 
finite temperature and important nonlocal contributions can be obtained, work for 
a thermodynamic system of massive fields as well. Besides, this is the most studied 
field model, so let us investigate an ensemble of mult i- component scalar massive fields 
satisfying equation, 

"(□ - -R - m 2 )l + P(<f))} up = 0. (57) 

LV g / J 

Because the mass term can be factorized out of the heat kernel, one can still use 
massless heat kernel (ET1) to derive the free energy, 

pP = _i r^ e sm^ K ^ s) (58) 

2(3 Jo s 



As usual, we subtract n = mode from the image sum, ([18]). Let us first treat 
local terms. The result in terms of the modified Bessel functions reads, 

+ ^A- 1 (m/3n)P + 0[» 2 ]}. (59) 

So far this expression is valid at any nonzero temperature. However, we are able to 
proceed and obtain explicit formulae only in high temperature limit. Simple expan- 
sions of the Bessel functions at (3 — > with the subsequent n-sum evaluation produces 
known local contributions KR EO]. The total result for free energy of massive fields 



at high temperature looks like: 

FL-FZ = -/^, 1/2 tr{^i + ^(P— 2 i) 
+ ^\R lj l?(-A)R i n + R^(-A)R1 



32tt 2 

+ P 7 f(-A) J R + P 7 f(-A)p] + 0[3? 3 ]} , (3 -> 0, (60) 
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The computational procedure for second order terms is performed after Poisson 
resummation (]5T|). Applying Eq. ([52]) to the basic form factor of nonlocal free energy 
for massive fields , 



7 



r°° ds 


2 




3 (o,e-^) - 1 


lo s 





e -sm J(_ S A) 



(61) 



(vacuum contribution subtracted in ([H]) is dealt with at the end of the present sec- 
tion), and using the integral 

/' 

Jo 



da(m 2 - a(l - a)A) ( ~ 1/2) 



2 f ^A 
arctan 



2m 



(62) 



we get 



7 (/V-A) 



4vr /V=A 
arctan 



871 



+ 



2m 

oo 

arctan 



k=l 



— dn arctan 
Jo 



Am 2 + 16rr 2 k 2 /(3 2 



(63) 



/Am 2 + 167t 2 k 2 /(3 2 i 

This equation is valid at arbitrary finite temperature, therefore, free energy of massive 



fields is nonlocal at any temperature. The first term of Eq. (|63|) came from k = mode 
of the sum, and it is nothing but the leading term of high temperature expansion, 
(3 — > 0. The difference of two divergent terms in the square brackets is finite, however, 
we are unable to give the result in a closed form. Thus, we restrict consideration to 
leading terms of high temperature expansion and understand the basic thermal form 
factor as, 

4vr / v /3 A* 



V(-A) 



arctanf ) + Of 

v 2m 1 



(3 -> 0. 



(64) 



The main nonlocality is contained in the leading term (pi]) . Subleading terms com- 
bined with vacuum contributions are not important at high temperatures. The full 



table of form factors 71(5), i = 1, ...4, in terms of q 



-A 



reads 



7i(9) = 75— 
pm 



— Q -?+(! + Aq 2 Y arctan — 
12 K 2q 



(65) 
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72(g) 

73(g) 
74(g) 



13 11 

— q - 3g 3 - (y + 2 8g 2 + ^8q 4 ) an-tan 



nq 

irq 

2(3m 

^arctanl-l. 
pm v 2g 



(£) 



2g + (1 - 4g 2 ) arctan(— ) 



(66) 
(67) 
(68) 



For practical purposes of physical applications we need to know spectral form 
representations for (|65|) — (j68|) . A spectral form for the basic form factor ([54]) is obvious, 

1 



/(-A) = f 



dm ■ 



(69) 



4m 2 - A ' 

Its massless limit immediately gives (0). Spectral forms for form factors with sub- 
tractions are obtained similarly (see appendix [B]). Then, all form factors (|BT)|)-(|68|) 
admit the form, 

1 



(70) 



where mass spectral weights are given in the table, 



pi(m 2 ) 



2m 2 m 4 



m 2 



rrr 



p 2 (m 2 ) 



p 3 (m 2 ) 



p±(fh 2 ) 



1 /25 14 m 2 



m 



32 V 9 
1 m 2 
3 + 



3 m? 



+ — 



rrr 



~ 9 ) 



(71) 

(72) 

(73) 
(74) 



Now we need to complete our derivation with the regularized free energy at zero 
temperature F°°. Nonlocal effective action for massive fields in an arbitrary space- 
time dimension has been calculated first by Avramidi E|. His approach is a direct 
summation of derivatives in a massive field theory, but we can make use of the mass- 
less heat kernel ([H]) obtained with the covariant perturbation theory and arrive at 
the same result. We compute zeta function regularized effective action according to 
the equation, 



W c 



ld_ 
2~d~e 



ipt poo Ac 



T(e) Jo s 



,1-e 



TrK(s) 



(75) 



<E = 
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Then, we get the following result for zero temperature free energy (the specific form 



of the effective action in Ref. [|2] in four dimensions), 
' - / <V:r <) 1 -tr \ 



32^ 



+ m 2 (in (^) - l) P + [RiM-&)R ij + Rl2(~A)R 
+ Ay 3 (-A)i2 + Ay 4 (-A)P] + 0[K 3 ]} , 
where form factors 7, are given in terms of dimensionless argument q = 



'—A 



(76) 
by the 



following expressions, 
7i(?) 
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-2 l + 4g 2 ' 



46 56 
15 + "3 

arctanh 



■<f + 32g 4 



(77) 



72(g) 



180 



97 
30 



- 17g 2 - 12g 4 



+ 2Jl + V(l + 8g 2 + 6g 4 )arctanh 
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73(g) 



74(g) 



- - 4g 2 + 4g 2 yl + 4g 2 arctanh 
3 

^ m (^") + 1 - \A + 4g 2 arctanh ( 



.VI + 4g 2 , 

1 



(78) 
(79) 
(80) 



This effective action may look more similar to Eq. (O) if the inverse hyperbolic 
tangents in functions 7 j are expressed in terms of logarithms. We have to remark 
here that form factor 73 is different from the others. Similarly to that of massless 
fields it does not depend on the regularization parameter \i. However, Eq. ([79]) is 
nonlocal in contrast to local fl4"6]). Of course, in the zero mass limit Eqs. (|76])-(PUD 
turn to Eqs. (H)-©. 

Finally, it should be noted that in order to use Eq. (|76|) for specific physical models, 
the convenient way to work with form factors (f77[)— fl8"0"|) is to treat them in a spectral 



form representation ||39|| . Then, the following mass spectrum integral is to be used, 
4 1 



V4m 2 - A 



arctanh 



V4m 2 - A / 



dm 



m 2 — m 



2 4m 2 - A ' 
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5 Free energy of massless spinor fields 

In this section we consider the massless Dirac spinors ip in a Euclidean ultrastatic 
spacetime at finite temperature. The massless covariant Dirac equation is taken as 

= 0, (82) 



where the standard notation )7 = 7 M V M is used (see |2T| for general definitions). The 
method of calculation of the effective action for spin-1/2 fields, W(i/ 2 ) , is similar to 
the one for spin-0 fields. The main difference is that fermions are antiperiodic in the 
Euclidean time and, therefore, they satisfy boundary conditions (fj) with the minus 
sign. The local form of the one-loop effective action W(i/2) = — Trln was studied 
first in Ref. H. It is defined in terms of the heat kernel (or zeta function) of operator 
([82]), however, following De Witt's idea we consider the squared operator y 2 , thus, 

Wfi/ 2 ) = Ki/ 2 ) = Tr^f 1/2) ( S ), (83) 

with the heat kernel Kn/ 2 )(s) corresponding to the squared Dirac operator. It can be 
shown || that the heat kernel (Green function) of the operator ft 2 is equivalent to 
the spinor heat kernel which is a solution of the equation, 



□ - -R{x)\ K (1/2) (s\x; x') = 15{s)S(x; x'). (84) 



One can represent the heat kernel (|8~4] ) at finite temperature in a form of the image 
sum 



43, 29 



£f 1/2) (s | t,x;t',x') = £ (~l) n K (1/2) (s | T,x;r' + f3n,x') (85) 



(cf. Eq . (|i2|) ) . Because time dependence of the heat kernel in ultrastatic spacetimes 
factorizes out (|T3"D, the trace of the heat kernel can be written in the form 

Tr#f 1/2) ( S ) = 2 ( O ,e-^)^p J d s xtrK$ 2) (s\x,x), (86) 

where 9 2 is the Jacobi theta function. It is convenient to use the ^-functions identity 

2 (O, e~ z ) = 3 (O, e" 2 / 4 ) - 9 3 (0, e~ z ) (87) 
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and express the heat kernel trace at finite temperature ( |86| ) in terms of the Jacobi 
functions 63. Thanks to this fact [29. 43] we do not have to repeat all calculations 
and can get the free energy of Dirac spinors at finite temperature using mathematical 
derivations of section |[ Spinor thermal form factors 7(1/2) are obtained then by a 
simple combination of scalar form factors 7 

7(1/2) (P\ 



-A) = - 



(88) 



To form the operator of Eq. (^4j) the potential term should be taken P = —j^Rl. 
The commutator curvature in not zero when covariant derivatives act on spinors, 
but we need to know only that tilZ^lZ^ = —^R^apR^^trl, where the squared 
Riemann tensor must be expressed via Ricci tensor and scalar curvature with help of 



the Gauss-Bonnet identity [44]. All other matrix structures are reduced to trl = 4. 
Taking into account these properties the free energy of massless spinors reads (from 
now on we omit spinor indices (1/2)) 

7tt 2 1 



/dW /2 { T 



80/? 4 144/? : 



:R 



+ 



% 7 f (-A)iF' + Ry$(-A)R] + 0[R 3 }} 



where the thermal spinor form factors 

2n 



_(3^At sh(2vrt/(/3v^A)) 
with the trigonometric polynomials 



1 

T 2 



9i{t) 
92® 



1 ^cos(t) ^sin(t)^ ^ 3 ^cos(t) sin(t) 



4V t 
1 

16 



t 2 



2 V *3 



t A 



sin(t) cos(i) rcos(t) sin(t) 



t 3 



(89) 

(90) 

(91) 
(92) 



t 2 t 

Note, that the only difference of (|90|) from ( |38| ) is the hyperbolic sinus instead of the 
hyperbolic tangent. 

This result is to be combined with the regularized contribution F™ n of the image 
sum (|85|) . The zeta regularized effective action ( j33f) at zero temperature takes the 
form, 

1 



57T 



d 3 x g 1/2 tr{R ij j 1 {-A)R ij + R l2 (-A)R + 0[K 3 ]}. 
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(93) 



where zero temperature spinor form factors 7i(— A) 



1, ...2, are 



7i(-A) 
7 2 (-A) 



1 

40 



In 



12 

y 



i 

120 



77 
30 



(94) 
(95) 



The main result for the finite temperature free energy Ff en is the sum of Eqs. (jS9|) 
and (|93|). It is an essentially nonlocal functional. On the other hand it is well known 
that only local terms survive in high temperature limit [53]. To check this property 
we consider the high temperature limit of our result. 

The calculation of (3 —>■ asymptotic of the thermal form factors (0) are anal- 
ogous to similar derivations in section |3j. A different hyperbolic function appearing 
in thermal form factors, namely, the hyperbolic sinus, results in the absence of linear 
in temperature nonlocal terms in high temperature expansion, for example, the basic 
form factor (g%(t) = — 2sin(t) in fl9H|)) reads 



7(/V- A ) 



-2 
+ 



C + In (J- 

7 « 3 V< 



TT 



-A) 



1 



31C(5 V(-A) 2 + 0[/5 6 ], (3 



0. 



(96) 



24tt 2 v ' 640?r 4 
Uniting the vacuum free energy (|9~3" ) and the high temperature expansion of (jS9|) into 
one expression, we observe that the logarithmic nonlocality also cancels in the sum. 
The resulting T = 1/(3 — > oo expansion looks like 



dW /2 
1 



7tt 2 



180/? 4 144/3 2 
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167T 
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— RR RaR^ 

30 10 3 



+ (3 



128 TT 4 
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1024 7T 6 
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RAR - —/>',; A /)"•' 
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-RA R 



—R ii A 2 R ij 
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+ 0[i? 3 ] + 0[/3 6 ]}, /3^0. 



(97) 



Combinations of quadratic in curvatures terms in the square brackets are the func- 
tional traces of Schwinger-DeWitt coefficient a 2 , a 3 , and a 4 after substituting the 
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Riemann curvature with the Ricci tensor and scalar curvature ||, [10], [38]. Expres- 
sion 



coincides with the massless limit of a similar expression in Ref. [29], which 



corrects apparent misprints in higher /5-order terms of Ref. [43 



Massive spinor fields can be treated in a similar way. At finite temperature the 
corresponding effective action is nonlocal. We do not present the result here because 
of its complexity. In the physically interesting limit of high temperatures the effective 



action of fermi fields becomes local and can be found in Ref. 29 1 . A small correction 



should be made to the result of |29j because the effective actions corresponding to the 
first order operator (J7 + m) and the squared one differ by local terms (see Eq. (4.2) 



of Ref. 45 



6 Conclusions 

Free energy of quantum fields in ultrastatic (optical) spaces has been a subject of 

]. When effects of gravitational 



study in a large number of papers []lj], [J0[ [43], [46], 
fields are negligeable, the metric under consideration is flat, and the spacetime is 
automatically ultrastatic. But even this simple case has not been studied sufficiently, 
as most works are concerned with finding the effective action either on constant 
background fields or (and) at very high temperature j|, |5|, [40]. However, rapidly 
oscillating background fields generate nonlocal terms in the effective action which 



contribute to vacuum polarization effects [g2|, |Tj, [48], [4JJ. Needless to say, that an 
interesting and important high temperature limit of the effective action is still only 
an asymptotic, and knowledge of the effective action behavior at arbitrary finite 
temperature is necessary. 

We have obtained nonlocal structures of the one-loop Euclidean effective action 
and free energy for thermal fields in asymptotically flat ultrastatic curved spacetimes. 
For nonconformal massless scalar fields this expression has been found at arbitrary 
finite temperatures. Explicit formulae for the high temperature limit of this general 
expression have been obtained. For massive scalar fields the free energy is derived 
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in the high temperature limit. With help of these results we calculated also the free 
energy of massless spinor fields. 

The calculated one-loop Euclidean effective action is known to be a generating 
functional of one-particle irreducible Green functions || |5(J. Therefore, variations 
of the nonlocal effective action over background fields generate the Green functions 



51] and the energy-momentum tensor [48|, while variations of the free energy over 



thermodynamical variables, such as temperature, provides one with thermodynamical 
potentials, entropy, etc., [[[£], pi]] . 

We would like to emphasize again that the technique we adopted, namely, the 
nonlocal covariant perturbation theory fl6|| , is a limit opposite to the effective po- 
tential method |5], [8] and the derivative expansions H 47, 52] suitable for constant 



or slowly changing background fields. Therefore, to make comparisons one has to 
expand our nonlocal results in powers of the derivatives and compare them with the 
corresponding terms in the known derivative expansions expanded, in turn, in powers 
of field strengths. Our results are in agreement with all known expansions of this 
kind. 

It is important to stress the value of the nonlocal effective action approach. As 
has been shown in section [3] the effective action for massless fields is naturally infrared 
finite, which demonstrates that infrared infinities are absent both at zero and finite 
temperatures. As for the applications of the covariant perturbation theory to other 
models, one can refer to Refs. |53], [54]], where spinor quantum electrodynamics in 
flat spacetime and scalar electrodynamics with gravity in the context of conformal 
anomaly have been studied. 

We used in this paper the imaginary time formalism, which means that the notion 
of the global periodic time is introduced. This greatly reduces the class of physical 
systems and spacetimes that can be considered. Alternatively, there exists real time 



formalism |55|] which can give a covariant form of the partition function |56[ and free 
energy. To work out nonlocal structures of the effective action using the perturbation 
theory, which is applicable to nonequilibrium systems and covariant from the outset, 
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one should derive it from the scratch. It would be an extremely interesting project 
to do. 

Obvious next step is to apply conformal transformation technique |28, |29[ to the 



obtained results and find the nonlocal free energy on a more general static gravita- 
tional backgrounds. This problem is facilitated by the fact that conformal transfor- 



mations of the nonlocal effective action are studied in detail in Refs. [p7| . In order 
to study particle creation by external gravitational fields the third curvature order is 
required [[23|, [49| . Our intention here was to develop a general method for calculations 
of nonlocal free energy, which is the most transparent while working in the second 
order in curvatures. Derivation of the next perturbative order, if needed, will not 
pose serious technical problems. 



Effects of space boundaries are also interesting and important |19, f|6] and they 



deserve investigation, however, to do so by means of curvature expansion a major 
revision of covariant perturbation theory is required. 

The results of this paper can be generalized to the very important case of back- 
ground gauge fields. Thus requires more caution in dealing with temporal components 
of gauge fields and leads to the appearance of other noncovariant terms in the effec- 
tive action [ESI. It is not surprising, since at finite temperatures Lorenz invariance is 



broken. But the general structure of nonlocalities in the effective action is still the 
same. 
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A Massless form factors with subtractions 



Complex form factors entering the trace of the heat kernel (^T|) are to be treated 
similarly to the basic one, as is described in section || The heat kernel form factor 
with one subtraction generates a thermal form factor in the effective action by the 
following computational procedure, 



f(sA) - 1 
sA 



(9* 



Using integral representation, 

~ 1 = /doWl-ai) / dase^ 1 -" 1 ^, 
sA Jo Jo 



(99) 



and taking the integral over proper time, we reduce it to the form, 

771(2) — 4 / d«i «i(l — ai) / d«2 Kq{hz ^/a2 J a\{l — cti)), (100) 
Jo Jo v 

where z = /3a/— A is a new dimensionless variable. This is a table integral with 

respect to «2, 

£dyy 2 K (by) = -^i(fo) + ~ (101) 



which, after introducing a new variable, y = 2Ja(l — a), give us 



n=i 



If 1 nzy, 
dy^===Ki( 



Jo a/1 — 2/ 2 



n 2 z 2 



(102) 



Applying the Bessel function relationship, 



^(2/)= - —K (y), 
dy 

to ii'i and using ( ^5|) for the resulting integral, we get, 



= - 4 E 



n=l 



dt 



sin(t) 



(103) 



'o (t 2 + n 2 z 2 / A) 2 n 2 z 2 
After integration by parts, to reduce the power of the denominator, rji admits the 
form similar to Eq. (|35|). 

/ sin(t) cos(t)\ 1 



— roa 
n=l J ° 



\ t2 
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(104) 
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With use of (|36|) it gives the final answer, 

"2tt 



I sin(t) cos(t) 



t 2 
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1 



1 



*t th(2irt/z) t 2 



Let us treat now the form factor with two subtractions, 



r /2 (/?^A) = ^ O °^[^(0,e-^)-l 



f(sA) -1-lsA 
(sA) 2 



(105) 



(106) 



In the representation, 

;(_ sA ) -1-hA 



sA 



d^a^il-a^ 2 / da 2 a 2 / da 3 e ai{1 ~ ai)a2a3sA , (107) 
jo jo 

it can be integrated first over the proper time. Integration over a 3 is exactly Eq. ( |101|) , 
and integral over a 2 -parameter is 

dyy 2 K 1 (by) = ~K 2 (b) + ^. (108) 

Then, 7] 2 looks in variables y and z as 



4E 
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n 2 z 2 Jo 



1 , V 3 ,nzy 1 1 

dy 7= „ ir 2 (-^-) + 



1 



3 n 2 z 2 



(109) 



It is possible to reduce (|109|) to 7 and rji types of integrals employing the relation, 



K 2 (y) = -K 1 (y) + K (y) J 

y 



(110) 
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12- 



n 4 z 4 



We already know how to deal with t-integral, what is new here that is a sum, 
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+ 7777 



7T 
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£ n 2 ^ 2 (t 2 + n 2 ^ 2 /4) 6z 2 t 2 8t 4 4^t 3 th(2vrt/z) ' 
which brings r^ 2 to the form, 
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dt I sin(t) + 3 
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And the final result reads, 

1 f°° 
= " 2 X dt 
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B Massive form factors with subtractions 



During the computation of remaining two form factors, with one subtraction 



2-y/? f°° ds 



f(sA) 



(3 Jo s 3 / 2 
and with two subtractions 

f°° ds 



sA 



f(sA) 



0-0, 



(115) 



/3-0, 



(116) 



(3 Jo s 3 / 2 (sA) 2 
we get rid of the image sum from the outset and keep only k = term. 

To perform the proper time integration we need two singular integrals that are 
regularized by cutoff at the lower limit, 



, f°° ds _ 



sb 



,5/2 



-2Vtt6 + 

4 ,3/2 i- 2 2b 



;ii7) 

'118) 



where b > 0. The form factors (|115|) and ( |116|) take then the form, 



Vi(P^A) = j£{j Q daA 1 {b{a))-A 1 {m a )), 



V2((3^A) = ^( £daA 2 (b(a)) - A 2 (m 2 ) - ^(m 2 )), 



(119) 
(120) 



where 6(a) 



m 



Aa(l — a). Now we make use of regularized integrals A\ and 



Ai- Functions r\\ and 772 are regular, therefore, all intermediate singularities in these 
equations mutually cancel. Integrations over a-parameter do not pose a problem and 
they are similar to the integral (|5^). The only difference is that higher powers of 
integrand b[a) bring extra factors at the arctangent function: 
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Spectral forms for r)i and 772 are readily obtained with help of basic spectral integral 



(|69|) applied to arctangent functions in ( 121 )— ( 122 ). To transform operator factors to 
spectral weights we only need to know the identity, 
11 111 



A 4m 2 - A 4m 2 A Am 2 Am 2 - A ' 
The spectral representations read, 



(123) 



,„ l — 7T r°° , / 2m 2 m 4 \ 1 , 

r?2 /V-A = — / dm(l ^ + ^)^ r 125 

' V V 7 4/5 im \ m 2 m 4 ^4m 2 -A v ; 
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